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ANSTHACT 


In  order  to  correlate  data  obtained  in  tests  at  the  Ordnance 
Aerophysic*  Laboratory  on  Von  Karman  radomes  with  theoretical  analy¬ 
ses  of  thermal  stresses  in  ogive  radomes,  the  temperature  and  flow 
conditions  around  the  none  of  a  blunt  radome  were  investigated  and  the 
results  are  reported  herein. 

The  blunt  noae  causes  a  detached  shock  wave  to  form  with  a 
transonic  flow  region  developed  between  the  wave  front  and  the  nose. 
Approximate  solutions  for  the  velocity  components  in  the  highly  nonlinear 
behavior  are  obtained,  and  detachment  distances  compared  with  published 
experimental  data.  Subsonic  and  supersonic  compressible  flows  are 
studied,  and  the  incompressible  case  is  shown  for  comparison. 

Two  numerical  examples  are  given  to  illustrate  the  evaluation  of 
initially  unknown  exponents  and  the  computation  of  constants  in  the  ve¬ 
locity  equations. 
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INTRODUCTION 

The  OAL*  aerodynamic  heating  tests  (Reference  (1)'  were  conducted  on  Von  Karas-n 
radomes.  Theoretical  analyses  have  been  developed  for  ogive  radomes  (Reference  '2V) 
for  the  prodiction  of  thermal  stresses.  Temperature-dependent  materiel  properties 
were  Investigated  In  Reference  (3)  and  radome  characteristics  In  Reference  (4). 

Based  on  these  investigations,  it  was  concluded  that  the  principal  factor  in 
expression  of  the  different  equations  and  correlations  between  test  and  theory 
is  rhe  coordinate  ar.^lo  (»),  which  is  defined  in  Figure  1. 

The  authors  of  Reference  (3)  have  pointed  out  that  there  are  no  theoretical 
solutions  for  flow  conditions  behind  a  detached  shock  wave  and  along  the  surface 
of  a  blunt  body,  but  they  do  present  experimental  data  that  is  plotted  In  Figure 
2  for  a  sphere.  The  test  measurements  were  obtained  for  the  ratio  (hQ/Ra)  of  the 
detachment  distance  divided  by  the  spherical  radius. 

Since  1957,  a  great  many  investigations  have  been  successfully  conducted  on 
high-speed  flow  problems  (References  (6)  to  (16),  inclusive).  Russian  scientists. 
References  (6)  and  (7),  set  up  a  system  of  integral  equations  that  they  solved 
numerically  by  extensive  computer  programs.  These  procedures  were  adopted  and 
extended  in  the  USA  between  1958  and  1963,  as  discussed  in  References  (8)  through 
(12)**.  In  the  field  of  physical  chemistry,  Lighthill  (13)  and  Whitman  (14) 
developed  formulas  for  shock  detachment  distances  based  on  gas  ionization  and 
dissociation  phenomena.  Later  (1964),  experimental  techniques  for  measuring 
detachment  distances  were  described  in  Reference  (15)  and  numerical  methods  for 
the  calculation  of  flow  properties  were  again  reported  in  1965  (Reference  (16);. 

The  present  study  grew  out  of  the  need  for  analyical  functions  in  thermal  - 
stress  analysis.  The  investigations  reviewed  in  the  preceding  paragraph  lead  to 
valuable  point-to-point  numerical  d<<ta  on  flow  characteristics;  and  the  functional 
forms  gained  herein,  usually  by  the  insertion  of  empirical  coef f icients ,  furnish 
mathematical  expressions  that  are  applicable  along  a  radome  surface  in  terms  of 
the  radome  coordinates,  which  govern  the  tbermal-st  ess  behavior. 

The  relationships  are  especially  necessary  at  and  near  the  radome  nose  when 
the  maximum  temperature  occurs,  and  the  thermal  stresses  tend  to  increase  owing  to 


*  Ordnance  Aerophysics  Laboratory,  General  Dynamics. 

**  An  excellent  bibliography  is  given  in  Reference  (12). 
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the  small  curvature  radii  that  prevail  in  this  region.  The  componental  velocity 
equations  are  presented  and  discussed  in  the  ensuing  text,  and  they  are  derived 
in  the  appendices.  Two  numerical  examples  are  worked  out  and  the  results  tabulated 
in  the  final  appendix. 


NOMENCLATURE 


The  following  notations  are  used  in  the  text  and  appendices  of  this  report: 


A,B,C,D 

F 

J,K 

M 

P 

R 

T 

U 

V 

X,Y,Z 


Constants 

Function 

Constants 

Mach  number 

Function 

Radius,  inches 

Temperature,  °F  or  °R 

Function 

Velocity,  fps 

Cartesian  reference  axes 


a,b 


c 

h 

k,m,n 

P 

q 

r 

x 

a 

0 

Y 

X 

P 

* 


Cons tants 
Sonic  speed,  fps 

Displacement  of  a  detached  shock  wave,  inches 
Exponents 

Pressure,  psi  or  psf 

Exponent 

Radius,  inches 

Auxiliary  variable  (x  ■  R/Rfl) 

Angular  coefficient,  radians 
Shockwave  angle,  degrees 

Ratio  of  specific  heats,  cp/cv,  (y  *  1.4  for  air) 

Constant 

Density,  scf 

Coordinate  angle,  degrees  or  radians 
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The  following  symbols  appear  as  subscripts: 
a:  Values  at  surface  "a" 

n:  Normal 

o:  Free-stream  values 

p:  Pressure 

s:  Stagnation 

t:  Tangential 

v:  Volume 

w:  Shock-wave  values 

R:  In  radial  direction 

t|l :  In  i|f  direction 


DETACHED  SHOCK-WAVE  LOCATION 


Approximate  Formal as --Hie  detachment  distance  hQ  (Figure  1)  is  the  distance  along 

the  geometric  Z-axis  between  the  shock  wave  and  the  radome  nose.  It  is  contained 

in  X  as 
w 


X 


w 


(1) 


where  R  is  the  radial  distance  from  the  origin  to  the  wave,  and  R  is  the  radius 
w  a 

of  a  spherical  body.  Therefore,  is  a  dimensionless  ratio,  and  the  variable 
ratio  is 


X  -  R/Rfl  (2) 

with  R  being  the  radial  distance  to  any  point  in  the  flow  region,  and  R  is 
positioned  by  the  coordinate  angle  i|i  as  sketched  in  Figure  1. 

The  detachment  distance  (hQ)  is  contained  implicitly  in  the  relation 


5 


»H«  JOHN*  MOM11M  UNIVIMIt* 

AFPLIICO  I»HYHIC*  tAUlONAroM^ 


(mfl-nk)  f 

[  m  +•  n  J 


O) 


wherein  Che  exponents  (k,  m,  n)  ere  dependent  on  Che  free-etream  Mach  number 
(Mq),  the  radome  dime niton  R#,  and  the  air-property  ratio  y  -  1.4,  which  la  the 
specific-heat  ratio  (y  ■  c  /c  ).  The  densities  appearing  in  equation  (3)  occur 

r  ” 

at  the  radoma-noee  stagnation  point  (pw#)  And  immediately  behind  the  wave  (Pw0) 
on  the  Z-axia.  Methods  for  the  calculation  of  exponents  k,  m,  and  n  are  described 
in  the  appendices,  and  the  denaitlea  are  computed  from  the  well  known  formulae 
which  follow. 


<y+Du  mJ 

4VmJ-  2 (y- 1 ) 


(4) 


P 


wo 


(P) 


pq  (y+D  M* 
~  Ty-O’mJ' 


(5) 


The  densities  (4)  and  (5)  can  be  calculated  directly  when  Mq  ie  known  or  assigned 
and  with  y  -  1.4.  Calculations  involving  equation  (3)  are  normally  carried  out 
by  the  trial-and-error  method,  and  the  various  relationships  are  discussed  more 
fully  in  the  appendices. 

Experimental  Data--Teat  data  on  detachment  lengths  are  reported  in  Reference  (5) 
or  spheres,  circular  cylinders,  flat  plates,  and  flat  nosed  bodies  of  revolution. 
The  tests  were  conducted  by  the  NACA,  NOL,  CIT,  and  in  Japan.*  The  discussions 
in  the  present  report  pertain  to  the  sphere  only,  and  the  test  data  plotted  in 
Figure  2  were  obtained  from  Reference  (5)  p.  105. 

*  NACA  »  National  Advisory  Committee  for  Aeronautics. 

NOL  “  U.S.  Naval  Ordnance  Laboratory. 

CIT  =  California  Institute  of  Technology. 
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AIN  VKI.CK'ITY  ■QUATIONS 


IncomPigaaihla  ?Itw»»Tli«  equation*  for  tneompreaalhU  flow  arc  alven  for  purpoae* 
of  rowparlaon  (Appendix  A)  and  arai 


VH  *“Vo  eo"  *  I  1 


•  +V0  ain  f  1  -f 


*1 


(6) 


J 

.1 


where  X  ia  tha  variant  ratio  X  •  R/R#,  and  tha  velocity  that  ia  tangential  to 
the  apherical  aorface  la  shown  below, 


(Vf) 


R  -  R. 


1,5  V  Min  f 


(?) 


Subaonlc  Compraeeiblg  riov--The  velocity  components  for  subsonic  and  aonlc  (M^l) 
compressible  flow  are  taken  from  Appendix  R  aa  follows: 


V„  -  -  V  coe  i 
K  0 


[S 


+  A  coa%  sin  v(X-l) 


(8) 


VQa in  » 


{ 


u  3AR  cos8* 

iu  .  - L -  }  1  +  coe  o(X-l) 

w 


(9) 


*  "  t"  ■  01  -  X  -1  ' 

a  w 


R  U" 

w 


»u. 


(10) 


w 


where  sin  a  0>w"l)  ■  0,  cos  a  (X^-l)  ■  1,  and  U  is  a  power  function  of  R  expressed 
in  terms  of  the  variable  X  as  written  below. 


Ot, 
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The  ax  rename  (k,  m,  n)  employed  In  the  function  t)  dopand  on  the  free-atream 
Mach  r  mber  (M^),  tl\a  radome  dimenalona,  and  tha  ratio  y  of  tha  apecifle  haata 
of  air.  Svaluationa  of  tha  axponanta  ara  axplainad  In  Appendix  D,  and  numarlcal 
valuta  ara  obtainad  In  Appendix  I.  From  tha  lattar  appendix  with  Mo  •  1,  tha 
numarlcal  raaulta  a rat 

H0  »  1 i  m  •  0.1;  n«  1.6204;  k  ■  12.202  (12) 

and  the  valuclty  which  ia  tangential  to  tha  outar  aurfaca  of  tha  aphara  pictured 
in  Figure  1  ia 


I 


I 


V 


a 


<V, 


VQain  t  (1.28-0,4  coa*<r) 


(13) 


which  can  ba  comparad  with  tha  of  aquation  (7)  that  waa  darivad  for  tha 
conditiona  of  incompraaaibla  flow. 

Suoaraonic  Compraaaibla  Flow- -Whan  tha  fraa-atraam  flow  ia  auparaonic  (Mq>1),  tha 
componanta  of  velocity  batwaan  tha  ahock  wava  and  tha  aphara  (Figura  1)  ara 
davalopad  in  Appendix  C  in  tha  form  given  below. 


VR  "  -Voeoi  ♦ 


‘  jflj.  B  ain  o^  (X-l) 

U*“  +  toa*y 
w 


+  A  coaBl|(  ain  a(X-l) 


(14) 


Vosin  ir 


♦  R 


m  + 

U' 

w 


BR 


:#coa  (X-l)  3AR#co sB(t 


co3*T 


ot 


{ 


1  +  cos  o(X-l) 


(15) 


x  -  r-  ,  a 

Ra 


KR  IT  ' 
o  w 


VI’  ai  '  2"<Vl>’  A  ‘ 


(16) 


(Y+1)  M, 


TT. 


K  -  ■ - 

y+1 


(17) 
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The  trigonometric  function*  at  the  outer  boundary  (X  •  X^)  become  *ln  <v»  ■  '• 
co*  (*x  (X  •  !)  •  0,  ain  cr(X^*l)  »  0,  and  eoa  afX^-l)  “  -1;  and,  at  tha  Inner  boundary 
(X  ■  1),  the  elne*  van  inti  and  the  cosine*  equal  unity.  The  function  U  la  defined 
by  equation  (11). 

For  large  value*  of  M  ,  the  conatant  »  (equation  17)  can  be  neglected;  and 
numerical  reunite  from  Appendix  8  are:  (See  Appendices  C,  D,  and  B) 


*  6 ;  m  ■  0,927;  n  ■  1 .4 ;  k  "  1 

Va  -  Vo  sin  *  (1.129  -  0.0114  cos"*) 
TEMPERATURE  RELATIONS 


(18) 

09) 


Inasmuch  as  the  present  investigation  was  made  to  find  suitable  analytical 
functions  for  the  computation  of  thermal  stresses,  the  velocity  funct  ons  are 
employed  to  represent  temperature  distributions  over  the  radome  nose.  By  means 
of  the  air  energy  equation,  the  surface  temperatures  (Tr)  are  expressed  as 


T 

a 


T 

o 


T 


s 


(20) 


in  which  ■  (T)^^  *«o  *8  t*ie  sta8nation  temperature  at  the  radome  nose,  and 

Cp  is  the  constant -pressure  specific  heat  of  the  gas.  Consistent  with  y  ■  1 .4 
for  air,  we  have 


c  -  6,006  (fps)3 /°R 
so  equation  (20)  can  be  written 

V? 

Ta  "  Ts  12,012 


(21) 


(22) 
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and,  with  aquation  (15),  ona  can  find  V 


'  KU 

-v1"*  sHh 

a  w 


coa*f 


.  &A.-£9J.aJ, 

cr 


For  a  specific  radoma,  va  shall  lat 


and  wnta  aquation  (23)  in  the  following  manner. 

V#  -  Voain  f  (Ao  +  A}  sec8*  -  A^cos8*) 


(23) 


(24) 


(25) 


One  of  the  advantages  of  the  preceding  type  of  equation,  arising  with  radome  noses 
that  are  only  approximately  spherical,  is  that  the  A^  can  be  found  empirically 
rather  than  by  applying  equations  (24). 

By  placing  equation  (25)  into  (22),  we  find  that  can  be  expressed  as  a 
series  in  cos  * ;  i.e. , 


T  -  T 
a  s 


n»3 

an(cosi|f)2n 


(26) 


where  the  a  are  listed  in  Table  1  below, 
n 


Table  1  -  Series  Coefficients  a 

_ n- 


n 

a 

n 

-2 

< 

-1 

A,  <2Ao-A1  ) 

0 

Aq  -  2A1(A0+  Ag ) 

1 

2AS(A1-  Ap)  -  A? 

2 

Ag  (2Ap  +  Ag) 

3 

< 
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In  eho  manner  JuaC  described,  wo  can  obtain  integrable  functions  of  f  that  are 
applicable  to  the  theory  presented  in  Reference  (2)  for  radome  thermal  stresses' 
The  calculations  with  equations  (25)  or  (26)  must  in  all  cases  bu  restricted  to 
values  of  the  coordinate  angle  in  the  range  0  <  ^  ,  where  ^  is  the  Mach 

angle  computed  from 

cos  ^  -  1/MQ  (2  7) 

and,  in  general,  since  the  approximations  limit  the  analyses  to  small  angles, 
one  should  probably  not  extend  the  calculations  beyond  about  30°  even  when  ^ 
exceeds  this  value. 

DISCUSSION 

The  approximate  solutions  obtained  herein  have  been  checked  in  a  few 
instances  with  experimental  data  and  seem  to  ba  satisfactory  for  engineering 
calculations.  However,  the  limits  of  their  applicability  have  not  been  fully 
determined  at  this  time,  and  it  appears  that  they  should  not  be  employed  much 
beyond  ■  30°. 

For  calculations  with  angles  larger  than  ^  =  cos"1  (1/Mq)  when  this  angle 
is  less  than  30°,  the  analyst  would  need  to  omit  the  sec2\|t  term  that  occurs  in 
equations  (14),  (15),  (23),  (25),  and  (26)  and  fit  a  solution  of  the  kind  re¬ 
presented  by  equations  (8),  (9)  and  exemplified  by  (13)  and  (19)  to  the  interval 
M  i  30°. 

CONCLUDING  REMARKS 

It  is  tentatively  concluded  that  the  velocity  tangential  to  the  radome 
nose  for  which  analytical  expressions  are  derived  in  this  report,  can  be  used 
to  compute  surface  temperatures  that  are  suitable  for  the  prediction  oi  thermal 
stresses . 

Furthermore,  it  appears  that  empirical  coefficients  can  be  employed 
to  obtain  semiempirical  functrons  in  the  present  solutions  ir,  cases 
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APPENDICES 

Appendix  As  Incompressible  Flow-Ifte  equation#  of  fluid  flow  are  nonlinear,  and 
their  general  solution  is  unknown.  To  illustrate  the  procedure  for  obtaining 
approximate  solutions,  we  shall  first  consider  the  case  of  incompressible  flow, 
which  also  provides  a  limiting  case  for  useful  comparisons. 

The  boundary  conditions  (Figure  1)  are: 


CD 


(VR)  --V0=os*, 

R  ■  00 


(V  )  -  V  sin  *  (Al) 

”  R  *  09 


<v 


R 


cy 


R  -  R 

a 


V(\jr) 


(A2) 


and  for  axially  synmetric  steady  flow,  the  continuity  equation  is  as  follows. 


3(prRVR)  5  (p  rV^  ) 

§R  +  dt 


0 


(A3) 


In  polar  coordinates  (R,  t|f),  the  cross-sectional  radius  is  r  ■  Rsin  t|i ;  and  the 
density  is  nearly  constant  when  the  free-stream  velocity  Vq  is  small  enough  to 
assume  that  the  flow  is  incompressible.  Under  these  conditions,  one  knows  that 
equation  (A3)  can  be  solved  by  employing  a  power  function  of  R  multiplied  by  a 
sine  or  cosine  function  of  i) .  Therefore,  in  order  to  satisfy  the  boundary 
conditions  (Al  and  A2),  one  chooses 


VB  =  -V  cos  f  ~  •  -V  U1  cos  if 
R  o  T  dR  o  T 


(A4) 


where  U  *  U(R)  is  a  function  of  R  alone.  From  potential  theory,  we  have 


avR  a(RV%) 
W  "  &R 


(A5) 
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which  together  with  (A4)  leede  to 

V  U  sin  t 
V  ■  -ft  . . -  . 

t  R 

If  we  further  select  U  ee  «  power  function  in  the  form 

U  -  R  +  <Ra/n)(Ra/R)n  -  Rft(X  +  l/nXn) 


i) 

II 

:l 

<A6)  ( 

'  1 
;j 

(A7) 


i 


I 


where  the  veriable  ratio  R/Rfl  is  represented  by  X  »  R/Rfl  and  n  is  an  undetermined 
exponent,  the  first  derivative  is 


U'  -  (1  -  1/Xn+1)  (A8) 

and  the  boundary  requirements  (Al)  and  (A2)  are  fulfilled  by  equations  (A4), 

(A6) ,  (A7),  and  (A8). 

The  unknown  exponent  n  is  found  by  substituting  the  foregoing  relations  into 
the  continuity  equation  (A3)  with  the  constant  density  term  cancelled  from  this 
equation,  as  shown  below 

d_ 

at 

which  tells  us  that  n  -  2.  And  the  velocity  components  become 

VR  -  -Vqco *t  (1-1/X3),  -  +VQsinf  (1+1/2X3)  (A10) 

which  is  the  classical  solution  for  this  problem.  When  R  ■  3R#  and  4Rft,  the 
components  are 


~  R*  sint  cos*  <1-1/Xn+1) 


Ra  sinaf  (X+l/nXn) 

»  ’ 


(A9) 


(V_)  -  -0.963  V  cos  t» 

R  X-3 


(V.)  -  1.019  V  sin  * 

1  X-3 


(All) 
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(V  )  -  -0.984  V  cos*,  (V.)  -  +1.008  V  iin  f 

K  X-4  °  ’  X-4  0 


(A12) 


from  which  it  is  seen  that  the  values  required  by  (Al)  are  quickly  approached. 
Along  the  spherical  surface  (Figure  1)  where  R  •  Rfi,  the  velocities  are 


(V, )  »V  •  +1.5  V  sin  ♦ 
♦  X-l  8  ° 


(A13) 


and  it  is  observed  that  the  tangential  velocity  increases  from  zero  at  the 

stagnation  point  to  V  *  1.5  V  at  each  side  where  tjt  *  90°.  According  to 

s  o 

Reference  (5),  the  latter  velocity  approximately  equals  the  free-stream  velocity 
Vq  when  1/  =  90°  in  the  high-speed  compressible  flow  of  air. 

Appendix  B:  Subsonic  Flow--In  the  consideration  of  subsonic  flows  that  approach 
the  sonic  range  (Mq  -  1)  as  an  upper  limit,  we  employ  a  more  general  function  U 
than  the  one  given  by  equation  (A7).  It  is 


^a  X111  1 

U  -  ~ 

km  vn 

nX 


in  which  k,  m,  n  are  undetermined  exponents,  and  we  take  the  density  as 


p  -  P(R)  F  (*). 


where  P  is  a  function  of  R  and  F  is  a  function  of  l|r. 

At  a  finite  radius  R  ■  R^,  as  yet  undetermined,  the  boundary  conditions 


)  X  V 
'R-R  o’ 
w 


(Vr)  «  -Vocos  h ,  (V.)  -  +V  sin  l|f 

R-R  ’  R-R 

w  w 


and  along  the  surface  of  the  sperical  nose  of  the  radome ,  the  boundary  require¬ 
ments  remain  as  specified  by  equations  (A2). 
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The  energy  equation  states  that 


yfl  V" 

T  -  T  +  t“ —  ■  T  4-  — 2- 
s  2  c  o  2  c 

P  P 


(B4) 


where  Tg  is  the  stagnation  temperature.  The  preceding  equation  can  be  expressed 


as 


-a-  1  +  XY-.UM8 

T  2  c  T  2 

P 

and  with  isentropic  compression,  the  density  can  be  calculated  from 


(B5) 


jr  -  wv 

s 


_L  _i 

Y-l 


!  +  hzpt  Uy 


(B6) 


which  can  be  expanded  to  obtain 


i  -  £  ( (2v-nif 

n  9  8  48  + 


(B7) 


and,  since  the  velocity  components  contain  sine  and  cosine  terms,  one  sees  that 
these  terms  will  occur  in  expansion  (B7)  as  cos3*,  and  higher  orders  with  similar 
powers  of  the  sine.  However,  since  sin3*  -  1-cos3*,  we  can  confine  the  dis¬ 
cussion  to  the  cosine  terms.  For  simplicity,  we  shall  consider  the  first  two 
such  terms  and  take 


VR  ■  -  Vo 


V  sin  * 


3 


where  equation  (B9)  is  obtained  from  (B8)  by  means  of  (A5). 


(B8) 


(B9) 
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Some  abbreviation  is  gained  by  letting 


a  -  tt/(X  -1)  (BIO) 

w 

and  from  the  boundary  requirements  (B3)  and  (A2),  we  find  with  the  aid  of  (Bl), 
(B8),  (B9) ,  and  (BIO)  that 


<V 


x-x 


w 


-V  cos  i|/ , 


U 


w 


R  U  ' 
w  w 


<VR)  “  0, 
R  X-l 


Vosinij( 


a 

R  U7 
a  w 


6Acosa 

ot 


A 


RU" 
a  w 

a  U  ' 
w 


(Bl  1) 


(B12) 


(B13) 


wherein  the  last  expression  (for  the  constant  A)  is  secured  by  requiring  the  maximum 
radial  velocity  to  occur  at  R  ■  R^  and  tjr  *  0  according  to  the  physical  condition  at 
this  point. 

When  V  (equation  (B12))  reduces  to  V  at  t  *  90°,  the  surface  boundary 
ft  o 

condition  becomes 


rib  - 1  <B1*> 

a  w 

and  further  relationships  between  the  undetermined  exponents  are  derived  from 
the  continuity  equation  (A3)  as  shown  below.  In  order  to  separate  the  variables, 
we  employ  the  following  approximations: 


£1 
U  ' 
w 


+  A  cosst|/  sin  or(X-l) 


U'cos  'll  (1  4-  a  cos2i|f) 


U  ' 
w 


(1  +  a) 


-cos  l|f 


fS* 


(B13) 


(B16) 
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■  uv««  IMIH  *»•*,«*# 


Jl  .  ii «! 

Vo  R 


3AR  con"*  / 

IT  “  - O' -  (  1  *  co,w^X"^ 


wherein  a,  b  represent  mean  values  of  the  functional  ratios,  l.e., 


A  U  '  aln  a<X-l)l 
— jl 


U 


mean 


crAU* 

tt  ^  wjt 

2R.U. 


TT 


(817) 


b  ■ 


3ARU '  f  ) 

— +  cos  o(X-l)j 


*Ua 


mean 


(818) 


and  with  (B2),  (BIS),  and  (B16)  put  into  (A3),  there  is  obtained 


&_ 

&R 


PFK^U  'sin  2*  (1  +  a  cos8*) 

"  2  (l  +7) V 


k 


m  fin.8.*  SI  2  b-£9».an 

u* 

w 


(B19) 


in  which  the  variables  can  be  separated  with  each  group  equalling  a  constant  that 
is  hereafter  denoted  by 


i_  ,  d(PHfV)  _  2  (1^  a)  .  d[p  sinail)Cl  -  b  cosail))'l 

UP  dR  "  F  "aln  2^  (1  +  a  cos8*)  d* 


(B20) 


The  preceding  equation  leads  to  the  following  two  ordinary  differential  equations 
involving  the  undetermined  coefficient  X. 


d  CPtTU')  _  XU  dR 
PR3U'  “  R2 U  ' 


(B21) 


dfF  8in8il(  (1  -  b  cos^il;)!  x  X  costy  (1  t  a  cosai)  dill 
F  sins\l/  (1  -  b  cos2 ty)  (1  +  a)  sin\|f  (1  -  b  corf) 


(B22) 
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Integration  of  (121)  and  (122)  produces 


WK, 


/vm+n  , v  v  r 

JL  ln  (x  '1>xw 

kmn  /w  nvfn  .  s  vr 

L  0^  -1)  X 


(123) 


^atnwf  (1  -  bico*®7, 


X  ln(sintcacf  ) 
-  b)  L 


+  _^_|i±»  +  4j  ln  I  ‘  '  b.CO‘> 

2(1  +  a)  I  1-b  bj  (  1  -  b  eo»** 


(B24) 


wharain  ♦  >  o.  In  ordar  to  hava  the  logarithm*  remain  finite  at  X  ■  1  and  ♦  •  u, 
wa  muat  require  that 


X  ■  kmn,  X  ■  2(l-b) 


(B25) 


with  the  first  expression  deriving  from  equation  (B23)  and  the  second  from  (B24). 
The  results  for  P  and  F  than  are 


,  u  iiiri  i 

m  +  n  X 

*  ■  w*’0'*1"”0  *mK/\ 


(B26) 


F  -  F, 


1  -  b  coa8^ 
1  1-b  cosBf 


(B27) 


where  the  exponent  q  is  evaluated  from  equations  (B24)  and  (B25). 


q  -  2  + 


Rl+a) 


(B28  ) 


From  equation  (B2)  the  density  equals  PF,  and  P  and  F  are  found  from  the  free- 

W  1 

stream  density  at  ii  ■  0,  X  ■  Xy. 
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k-1 

■ 

(B29) 

To  summarise  the  conditions  tor  computing  the  unknowns,  w«  find  by  means  of  (Bl), 
(BID,,  and  (B23)  that 


p  -  PF  *  p. 


Jl 


1  -  b  cosi)| 


(VX> 


(m+l-nk) 


m  + 


m+n 


m  + 


X"*"  •  (130) 

and  for  (B29),  whan  X  •  1  and  1)  ■  o,  tha  danslty  ratio  is  as  shown  balow  for  the 
stagnation-point  and  frae-etream  dtnsitias. 


k-1 


(m+l-nk) 


m  +  nX 

_ JL 

m  +  n 


m+n 


By  maans  of  (Bl),  (B14) ,  (B30),  and  (B31),  thara  is  obtainsd 


and  one  also  finds  from  (B30)  that 


,  „  m+n 

m  +  n  X 
w 

m+n 


X 

X-n 


(B31) 


(B32) 


(B33) 


which  can  be  put  Into  equation  (B31)  for  further  simplification  of  the  latter. 
The  set  of  equations  now  can  be  listed  as 
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(a)  X  •'  kmn 


(b)  X  -  2(1 -b ) 

<<=> 

W  »./»„  -  X„<—k>  [X/a-n)]1-1 

V  m+2 

,  .  X  Po  *w 

<•>  x^*"7 - 


(B34) 


and  the  solution  of  this  set  of  equations  (B34)  is  taken  up  in  Appendix  D. 

Appendix  C:  Supersonic  Flow-Ahead  of  the  blunt  body  in  Figure  1  for  supersonic 
flow,  a  detached  shock  wave  occurs  that  provides  a  boundary  surface  for  the 
definition  of  R^.  Air  properties  behind  the  wave  front  are  related  to  the  free- 
stream  properties  in  the  usual  way  through  application  of  the  continuity,  momentum, 
and  energy  equations  and  the  universal  gas  law.  Briefly,  with  n  denoting  normal 
and  t  tangential  components,  these  relations  are: 


p  V  -  p  V 
o  on  w  wn 


p  +  p  V  2  *  p  +p  v  2 
ro  ro  on  1  w  w  wn 


V  -  V 
ot  wt 


c  T  +  ~ 
p  o  2 


C  T  -f  —  m  C  X 
p  w  2  P  s 


p  *p  RI 

g 
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V  «  V  sin  p,  V  *  V  cos  B 
on  c  ot  o  M 


ca  -  Y  P/P  ■  <Y-1)  c  T 


the  foregoing  equations  lead  to  the  following  well  known  relations. 


Po  ^wn 

pw  Von 


2[1  +  (y-1)M0’  /2] 

<Y  +  1>  >*"„ 


the  following  constants. 


3  a  V+1  '-X  J(  a  B  m  ■  a 

B1  2  ’  * 1  Ki  2  *  °  B^P  * 

Therefore  by  insertion  of  (C9)  in  (C8),  we  may  write 


(C6) 

B 

(C  7  ) 

0 

.  And 

0 

(C8) 

!  ’■  | 
Lj  ! 

,  1 

employ 

Li 

li 

(C9) 

n 1 

V  1  +  K.  (M  sin  p)s 

wn  *  N  o 


B  M  3  sin  p 
l  o 


■  B  esc  P  +  K  sin  p 


(CIO) 


and,  with  y  ■  1.4,  the  above  constants  are 

B:  «=  6/5,  -  1/5,  and  K  -  1/6, 


The  general  boundary  conditions  along  the  spherical  nose  surface  (R  *  R  ) 

fl 


<VR)  -  o,  (V.) 

R  -  R  ’ 


R  -  R 


Va  =  V0|f  ,Ra) 


(Cll) 


and,  immediately  behind  the  shock  wave  (R  ■  Rw),  the  boundary  requirements  are 
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<v 


R-R 


Vwnsin(g-^)  +  Vwtcos(e+t) 


(V.)  «  V  sin(B+|)  -  V  cos(g+t) 

*  R-R  Wt  Wn 

w 


(C12) 


but  conditions  (C12)  are  intractable  at  the  present  time  owing  to  the  difficulty 
imposed  by  the  variable  Ry  and  the  unknown  analytical  relationship  between  (3  and 
i|i .  Some  progress  can  be  made,  however,  by  restricting  the  analysis  to  limited 
ranges  of  the  coordinate  angle  l|r .  If  we  consider  the  case  wherein  the  flow  region 
has  no  large  departure  from  the  longitudinal  axis,  we  can  use  8  ~  90°-ijr;  and  the 
simplified  boundary  conditions  become 


<V 


R-R 


-  V  =  -  V  (B  sec  i|i  +  K  cos  i) l) 
wn  o  f  T 


w 


(C13) 


(V, 


4) 


+  V  *  +  V  sin 
wt  o 


R-R 


w 


(C14) 


For  these  latter  requirements,  one  can  use 


*  TT^2(Xw"1)» 


or  =  tt/(X  -1) 
w 


(C15) 


V_  =  -  V  cos  i|i 
R  o 


B  sin  Qfj  (X-l) 


777-  + - 3  : . '  +  A  cos2^  sin  ar(X-l) 

U  cos*** 

w 


(C16) 


V  sin  ill 

V  =  -2 - 1 

ill  R 


BR  cos  Of  (X-l)  3AR  cos2l|f  f  ~] 

KU  .  a  i  a  /t.  -  /•>  •>  \  1 

+  - H  .  -  1  - - 1  1  +  COS  Of  (X-  1 )  ( 


U  '  a  cos8 iji 
w  1 


a 


(C17) 


where  U  is  the  power  function  of  R  given  by  equation  (Bl),  and  the  density  will 
again  be  represented  by  PF  as  shown  in  (B2) .  Again,  as  in  Appendix  B,  we  shall 
use  mean  values  of  the  functional  ratios  in  order  to  integrate  the  continuity 
equation  (A3),  and  we  shall  find  A  in  the  same  manner  as  described  below  equation 
(B13).  Hence, 
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V  M  - 


V  KU 'cos  ♦ 
o 

U’ 

w 


(1  +  a^sec3^  +  a  cos3^) 


V  KU  sin  ill 

V  ~  ~ - - - 

ili  R  U ' 
w 


(1  +  ^sec3^  -  b  cos  ill) 


A  ■ 


KR  U" 
a  v 

oU  ' 
w 


_BUJ  sin  ax  (X-l) 

B_ 

[1  +  53i  ) 

KU' 

<v 

2K  ! 

1  +  R  U" 

mean 

a  a 

AU  '  sin  a (X-l) 
_ w _ 

KU’ 


a  A  U 


w 


mean 


2  KR  U  ” 
a  a 


BU 


f  *  R  cos  or.  (X-l) 
w  a  1 


atKU 


BR  U  ' 
a  w 

2a,  KU 


mean 


3AR  U  ' 
a  w 

a  KU 


1  +  cos 


t¥<x-i)j 


3AR  U  ’ 
a  w 

a  KU 


mean 


Integration  of  the  continuity  equation  is  completed  by  the  method 
viously  in  Appendix  B,  and  we  learn  that 


X  *  kmn, 


P  cos  * 
wo 


2(1  +  bx-  b) 
i  +  al+  a 

k-1 


ci+  1 

V  1  " 

qa 

m  + 

„  m+n 
nil 

w 

Cj  +cosst 

A*  - 

Ca-cos3lf 

m  + 

„m+n 

nX 

<V*> 


(nrt-1-nk) 


(C 18) 

(C19) 

(C20) 

(C21) 

(C22) 

(C23) 

(C24) 

explained  pre- 

(C25) 

(C26) 
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P  B  Ma 
o  l  o 


p  (Y+1)M® 
0  o 


1  +  i^M3  2  +  (Y-l)M3 


(C27) 


p  r(w* 

W° 

L  1  O 


(C28) 


where  B^  ,  K  are  defined  by  equations  (C9)  and  ^ ,  Cs  ,  qJ ,  and  qg  are  given  below. 


C  =  -  JL  +  Zl_  +  _L  c  -  1  ,  l  ^  \ 

1  2b  +  y4bs  +  b  ’  Cs  ~  2b  +  ^4b»  +  ~ 


(C29) 


„  1  +  °i  \  f1  +  a<1“Gl)l  1  -  Ca  x  fl  +  ad+Cg)] 

ql  =  ‘  C1+  C2  -  2b[  C1+  Cs  J*  *  "  2  +  C^Tc^  +  2b[”c"+  c2 - < 


(C30) 


When  the  foregoing  results  are  summarized  for  evaluation  of  the  unknown  exponents, 
we  obtain  a  set  of  equations  that  are  similar  to  (B36) . 


(a)  X  =  kmn 


(b)  X 


2(1  +  bx  -  b) 
1  +  ai+  a 


y_\  y,  m+n  X  +  mK  \  +  m 
'  '  w  X-nK  “  X^"7 


(d)  p  /p  *  X 

KW8  HWO  W 


(m+l-nk) 


^\/(\  -n) 


.  Pvo^V' 

'1_n  Pws^l^X-1) 


(C3 1 ) 


The  above  set  differs  from  the  previous  one  primarily  in  the  use  of  X  *  >,/K 
(where  X%  -  6X  when  y  -  1.4),  the  density  ratio,  and  the  additional  terms  in 
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equation  (e).  The  density  ratio  is  calculated  with  equation  (C28),  which,  when 


Y  ■  1.4,  gives  us 


5 

2  1  2 

36Mo 

5C7M*-  1) 


(C32) 


and  pwq  can  be  computed  with  equation  (C27)  when  required.  Further  discussion 
of  the  solution  of  set  (C31)  is  presented  in  the  next  appendix. 


Appendix  D;  Solution  of  Equations--Inspect5on  of  the  equation  set  (B36)  indicates 
that  two  solutions  are  m  ■  n  «  0  and  n  ■  -m.  These  two  solutions  are  trivial, 
however,  and  we  must  look  beyond  them  to  find  a  usable  function  U  (equation  Bl). 

In  doing  so,  one  can  regard  as  a  known  quantity  based  on  the  experimental  data 
reported  in  Reference  ($)  and  plotted  in  Figure  2.  The  first  solution  then 
obtained  proceeds  as  follows:  If  we  use  the  mean  value  of  equation  (B18)  together 
with  (B13)  for  A,  we  arrive  at 


A 


orXX 

w 


b 


<ypo(X-n)Xwm+2 


3ps(X+m) 

p  Q^X  m+3 
o  w 


(Dl) 

(D2) 


and  die  known  density  ratio  can  be  replaced  by  an  exponent  (m^ )  of  Xw  as  defined 
*jy  equation  (D3)  immediately  below. 


In  (ps/po) 
^  "  2  '  tn  X 

w 


In  this  way,  we  have  the  relation 


P 


o 


(D3) 


(D4> 


and  the  equations  (B36)  can  be  written  to  furnish  the  following  set  wherein  we 
let  of  Xw/6  equal  the  known  constant  C. 
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(a) 

X  “  kmn ,  C  ■  o?X  /6 
*  w 

(b) 

X  -  2«-b>  -  2  - 

(c) 

„  m+n  X  +  m 

X  =  . 

w  X  -  n 

r  n  k”  1 

(d) 

1  =  X  ^m+mi -1"nk) 
w 

-11 ' 
_ 1 

(e) 

X  .  m-taj 

X-n  w 

(D5 

-e)  into  (D5-d),  we 

obtain 

it  « 

1 

m  +  -  n 


(D5) 


(D6) 


and  from  (a)  and  (b),  one  finds 


£  *  X(l+C)  +  m  -  n  -  2C 


(D7) 


which  leads  to  a  calculable  value  of  X  when  (D6)  and  (D7)  are  solved  together. 

m,-  2 


X  -  2  + 


1  +  C 


(D8) 


The  coefficient  X  is  less  than  2  (since  is  smaller  then  2)  as  computed  from 

the  preceding  expression.  After  computing  the  above  coefficient,  n  and  k  are 
expressed  in  terms  of  m  by  means  of  (D5-a)  and  D6). 


X  (m  +  TDj  ) 

n  ■  — r— - 

X  +  m 


X  +  m 
m(m  +  ) 


(D9) 


By  using  (D9)  to  replace  n  in  (D5-e),  we  find 
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(X-m  )  O  "  m  (DIO) 

l  w  x  m 

w 

and  m  can  be  calculated  from  (DIO)  by  successive  trials,  because  all  quantities 
but  m  are  known  in  this  relationship.  Following  the  numerical  evaluation  of  m 
from  (DIO),  n  and  k  are  computed  with  equations  (D9).  These  values  of  the 
exponents  (k,  m,  n)  furnish  sufficient  information  for  determination  of  the  power 
function  U  of  equation  (Bl).  The  solution,  however,  is  approximate,  because  (D5«c) 
can  be  rearranged  to  give 

(X-n)X  n  -  £X.»t-S>  (Dll) 

w  x  m 

w 


and,  by  comparing  (DIO)  and  (Dll),  one  observes  that  the  latter  is  satisfied  only 
by  n  ■  ,  which  fails  to  satisfy  all  of  the  other  equations  in  set  (D5).  The 

approximation  thus  formed  meets  all  of  the  boundary  requirements  in  Appendix  B 
with  the  exception  of  (B14),  which  becomes  greater  than  unity,  lying  between  1 
and  about  1.5.  This  result  may,  nevertheless,  be  acceptable  in  many  cases,  because 
it  appears  likely  that  (B14)  is  too  severe.  We  are  unable  to  answer  this  question 
conclusively,  however, since  the  experimental  data  are  insufficient  for  definition  of 
the  correct  relation. 


It  is  necessary  to  consider,  also,  that  perhaps  all  of  the  boundary  conditions 
be  relaxed  to  some  extent  as  indicated  by  the  solution  in  Appendix  A  for  in¬ 
compressible  flow.  This  solution  shows  that  the  air  properties  differ  from  the 
free-stream  values  for  all  finite  values  of  X^.  Employing  this  behavior  as  a 
guide,  one  can  reduce  requirement  (Bll)  by  a  reasonable  percentage,  such  as  is 
exemplified  by  equations  (All)  and  (A12). 

If  we  say  that  (Bll)  is  J-percent  of  unity  (for  example,  one  might  use 
J  ■  0.9  for  a  ten  percent  reduction),  equations  (D5)  become 


(•)  X 
(b)  X 


kmn, 


Xj*  JX 


2(1  — b ) 


’  (Xa-fm)(X8-n) 


+ 


w 

j 
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(c) 

^  m+n 
w 

a ' 

“rr 

a 

(d) 

p,/po 

.X  (l 
w 

X 

(e) 

3 

.  a  Y  1 

X  -  n 

a  A 

w 

(rn+l-nk)  f 
m+n'a 


V<Vn> 


k-1 


(D12) 


where,  in  equation  (D12-e),  depends  on  the  choice  of  J  as  defined  by  the  next 
equation  (in  which  nij  is  given  by  (D3)). 


"a  =  +  iiV"'"  (D13) 

w 

Upon  solving  the  set  of  equations  (D12),  one  finds  that  X  is  no  longer  independent 
of  m  and  n  as  was  the  case  with  set  (D5)  where  the  result  is  given  by  (D8). 
Instead,  we  proceed  as  follows:  With  D  given  by 

D  ■ 1  +  frir  <D14> 

w 

we  find  from  (D12-d)  and  (D12-e)  that 


k  ■  D/(m  +  n^-  n)  (D15) 

and  by  means  of  (a),  (e),  and  (D15),  n  is  obtained  as  a  function  of  m. 


n  -  to,  +  m^l  -  (l-l/Xwm4l^)JD 
By  substituting  (D12-e)  into  (D12-c),  we  find 


(D16) 


(D17) 


which  together  with  (D16)  is  solved  by  successive  trials  to  obtain  m.  With  the 
value  of  m  thus  found,  n  is  calculated  from  (D16)  and  k  from  (D13).  And  these 
three  exponents  furnish  sufficient  information  for  the  definition  of  U. 
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Similarly,  tha  aquation  set  (C31)  can  ha  evaluated  by  taking  (C31-e)  In 


tha  form 


(D1H) 


where  m,  la  computad  from 


<V  2  - 


4n<pwa/pwo^n^*l/Blwl  > 


(019) 


and  tha  expression  for  k  that  la  almllar  to  aquation  (013)  la 


k  »  0^  /(m  +  n) 


(020) 


and  the  n-function  of  m,  almllar  to  (016),  becomes 


n  -  n^+  m  1  -  /K 

o.  « 


(D21 ) 


In  which  the  term  is  calculated  from  the  next  axpraaaion, 


Di«  1  - 


In  (1  -  i/B^ot) 


(D22) 


Using  the  trial-and-error  method,  we  then  calculate  m  for  the  following  relation¬ 
ship  in  conjunction  with  equation  (021). 


1  +  -  1  - 

n  x  m+n), 

w 


(D23) 


It  is  observed  that  the  preceding  relation  is  analogous  to  equation  (017);  and 
after  m  is  computed,  n  and  k  are  determined  by  means  of  equations  (021)  and  (D2U) , 
respectively . 
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Whim  it  in  desired  to  employ  ( h*  induction  coefficient  J,  th«  exponents  (m, 
o,  end  k)  are  computed  from  th*  fnl lowing  sequence  of  equations. 


n  -  n^+  m  ^1  -  (1  •  1  )JD^  /K 

x""*  .  ,  +#[i  --X-  " 

w  n  X 

»  "  ■ 


(024) 


(025) 


k“D,/(m  +  n\l»n) 


(D26) 


Exponent  m  in  calculated  by  trial  and  error  from  equation*  (025)  and  (024)  followed 
by  the  evaluation  of  k  from  (026).  The  terma  and  are  obtained  from 


"W  "*4 


(027) 


V  I\  + 


*"*w 


(028) 


where  and  D1  are  given  by  equatlona  (019)  and  (022),  respectively. 

Appendix  K:  Numerical  Bxemplea--If  we  consider  the  case  of  M0  ■  1  and  tale  ■  4 
with  the  solution  represented  by  equations  (D21)  to  (D17),  inclusive,  and  J  »  0.9, 
the  initially  calculable  terms  are: 


P m/pn  -  1  .577  (Eq.  032),  n*  -  1.671,4/0  (Eq.  03) 

9  O 


n^-  1.595,605  (Eq.  013),  D  »  0.924,135  (Eq.  D14) 


n  -  1.595,605  +  m  1  -  0.831,722  (1-0,109,514  x  4-m)J  (Eq.  D16) 
4n~”^  -  1  +  JJ  £l  -  (0.109 , 514) /4mj  (Eq.  017) 
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Uao  of  the  trial-and-error  proceaa  furniahaa  approximate  fulfillmant  of  equations 
(D16)  and  (D17)  with  m  ■  0.1:  and  from  (D16)  and  (D15),  n  and  k  ara  computed  and 
aummariaad  balow  with  m. 

m  ■  0.1:  n  ■  1.6204;  k  ■  12.282  (Rl) 


Tha  power  function  (Bl)  with  axponanta  (SI)  ia  ahovn  below  (wherein  X  * 


R/Rg). 


10  x 


12.282 


Tha  first  and  second  derivatives  with  respect  to  R  become 


. 


11.282 


x0.9  x2.6204 


R  U'  '  -  11.282  10  X0,1  +  ^ 

a  xi. 


„1.6204 


10.282 


x0.9  x2.6204 


+  10  xoa  4  3.61Z.*i?2 

1U  *  x1.6204 


11.282 


2.6204  0.9 

x3.6204  ’  x1.9 


and  from  equations  (E2,  83,  and  4),  the  following  numerical  results  were  com* 
puted. 

U  /R  -  3.255,587  x  10* l;  U  /R  -  9.171,341  x  1011  (E5) 

A  A  W  A 


u '  -  o 

a 


U'  -  2.542,788  x  10 
w  ’ 
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*  -  VQcos  i| l, 


0.902  V  sint 

o 


(E13) 


and  at  the  inner  boundary  along  the  surface  of  the  spherical  nose  where  X  =  1,  the 
components  are  as  follows: 
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(V.)  -  0,  (V4)  -  V,  -  V  s in  S  (1.280  -  0.4  cos3*) 

R  X-l  f  X-l  *  0 


(£14) 


The  results  indicate  that  V.  Increases  from  zero  *t  |  ■  0  to  V  ■  1.28  V  at 

a  a  o 

*  -  90°  with  the  variation  tabulated  in  the  immediately  following  Table  2. 


Table  2  -  Variation  of  Velocity  Tangential  to  Spherical  Nose  for  Mq  *  1 . 


♦ 

(D"g.) 

sin  * 

cos  * 

V  esc  * 

B 

V 

0 

VVn 

a  o 

(compressible) 

V  /V 
a  o 

( incompre  s  s ib le ) 

0 

0 

1 

0.88 

0 

0 

30 

0,5 

0.866 

0.98 

0.490 

0.750 

45 

0.707 

0.707 

1.08 

0.764 

1.061 

60 

0.866 

0.5 

1.18 

1.022 

1.299 

90 

1 

0 

1.28 

1.280 

1.5 

If 

one  next 

considers  the  case  for  large  Mach  numbers  (say  M  — ►  *) ,  the  value 

of  *w  * 

1.143  is 

taken  from  the 

data  in  Figure 

2  and  PWS/PWQ  - 

1.073  computed  with 

equation 

(C32) . 

Since  the  following  ratio  is 

small;  i.e., 

(1/aiBj.  ) 

•  5/2 Itt  ■ 

M  >6 

0 

0.0758 

(E15) 

and  will 

be  neglected,  we  shall 

let  k  »  1  and 

calculate  m  and  n 

from  equatior..-, 

(C31-a, 

c,  and  d) 

.  We  find  that 

n  *  m 


+  1  '  t  &U3)  ■  "  +  °-472'659 


(E16) 


(1  143)"^"  - 
U,i4:5;  n(6m  -  1) 


(E17) 
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where  (E16)  and  (E17)  are  solved  by  successive  trials  to  obtain  m  -  0.927  and 
n  ■  1.3997.  The  power  function  of  R  (equation  Bl)  wherein  X  ■  R/R  is  then 


(E18) 

(E19) 

(E20) 


U  *  1.793.035R  , 

U  = 

1.813.667R 

(E21) 

a  a’ 

w 

a 

U  '  =0 
a 

U  '  - 
w 

0.264,660 

(E22) 

RU"  -  2.327, 
a  a 

R  U'  ' 
a  w 

=  1.461 

(E23) 

The  ratios  that  appear  in  the  expressions  for  the  components  of  velocity  were 
also  evaluated;  e.g., 


and  the  velocity  components  are  written  as  follows. 


(E24) 

(E25) 


33 


THC  JOHNS  HOPKINS  UNIVfSSITV 

APPLIED  PHYSICS  LABORATORY 

SlCVIN  SPKtNtt  MAHVUAND 


c0,  * 


f_LLL 

L 1.588 


+  (0.0418)  cosaf  sin  7tt(X 


Vo8int|f 


U 


,  5O0  -  (0.00571)  Rcosl  \  1  +  cos  7tt(X-1) 

i.joo  a 


{■ 


i 


(E26) 


(E27) 


<V„) 


-  Vfi  -  Vosin  f  (1.129  -  0.0114  cosat) 


(E28) 


The  variation  of  V  along  the  surface  of  the  spherical  nose  shown  in  Figure  1  is 
& 

tabulated  below  in  Table  3. 


Table  3  -  Variation  of  Velocity  Tangential  to  Spherical  Nose  for  Mq>6 


♦ 

(Deg.) 

sin  if 

cos  if 

V  esc  if 
a  T 

V 

o 

V  /V 
a  o 

(compressible) 

V  /V 
a  o 

(incompressible) 

0 

0 

1 

1.118 

0 

0 

30 

0.5 

0.866 

1.120 

0.560 

0.750 

45 

0.707 

0.707 

1.123 

0.794 

1.061 

60 

0.866 

0.5 

1.126 

0.975 

1.299 

90 

1 

0 

1.129 

1.129 

1.5 

0 

0 

3 

3 

3 

3 

3 

J 

3 

3 

3 

3 
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